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I.  Integral  Equation  Status 


The  required  inner  products  in  the  Galerkin  solution  of  the  integral  equation,  using  a 
basis  of  exponential- weighted  Laguerre  polynomials,  are  ex'^rcssed  as  series  of  algebraic 
forms.  We  are  presently  evaluating  these  series  to  get  a  feel  for  the  invertibility  of  the 
integral  equation  (transformed  to  a  system  of  linear  equations),  prior  to  any  inclusion 
of  the  anticipated  asymptotic  behavior  from  the  Sommerfeld  half-space  problem.  If  this 
unmodified  system  of  the  first  kind  inverts  nicely,  then  its  inverse  can  be  used  to  express 
the  boundary  value  problem  as  a  system  of  the  second  kind;  resulting  in  a  whole  new  set 
of  opportunities  for  iterative  or  perturbation  series  approaches. 

There  are  several  important  issues  to  be  addressed  in  connection  with  the  far  (r  — >  oo) 
behavior: 


(1)  The  far  behavior 


oo,y 


0)  ~ 


k^a 

\/Swiki  x^!'^ 


is  not  the  proper  form  close  to  the  apex.  An  abrupt  jump  from  the  Laguerre  series 
to  this  oscillatory  function  is  likely  to  lead  to  numerical  instabilities,  which  suggests 
that  a  scaling  funtion  that  filters  out  its  effect  for  small  r  could  be  useful. 

(2)  Inclusion  of  this  far  behavior  is  complicated  in  context  of  the  even/odd  symmetry 
decomposition.  One  simplification,  based  on  the  same  line  of  physical  intuition 
used  to  insure  uniqueness  in  these  types  of  scattering  problems,  is  to  invoke  a 
radiation  condition  inside  the  wedge  material.  That  is,  assume  the  wedge  to  be 
slightly  lossy  so  that  the  effect  of  the  image  source  (in  the  even  or  odd  case)  is 
negligible  (in  comparison  with  the  original  line  source),  which  is  the  domain  of  the 
surface  integral  equation.  Admittedly,  this  is  only  a  conjecture  at  the  present,  and 
its  utility  remains  to  be  tested  via  experiments  (both  analytical  and  numerical). 
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II.  Conferences,  Meetings,  Propos.\ls,  and  Papers 

(1)  We  submitted  an  abstract  titled  “Single  Surface  Integral  Equation  for  Penetrable 
Wedge  Scattering”  to  the  Acoustical  Society  of  America  Meeting  to  be  held  in 
Cambridge,  MA  during  June  1994  (copy  attached). 

(2)  We  will  present  a  seminar  on  this  research  at  the  Electrical  Engineering  Department 
of  the  University  of  Mississippi  on  3  March  1994.  While  in  O.xford,  we  will  also 
meet  with  acoustics  researchers  from  the  National  Center  for  Physical  Acoustics. 

(3)  We  submitted  a  proposal  “Analytical  and  Asymptotic  Diffraction  Theory  for  Pen¬ 
etrable  Rough  Surfaces”  to  ONR  through  the  DEPSCoR  program  to  continue  and 
extend  this  acoustic  scattering  research. 

(1)  The  wedge  paper  [62],  which  is  a  result  of  last  year’s  grant  and  has  been  the  basis 
for  our  present  grant,  has  appeared  in  print  (copy  attached). 
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Single  Surface  Integral  Equation  for  Penetrable  Wedge  Scattering. 

Anthony  M.  J.  Davis  (Mathematics  Department,  University  of  Alabama,  Tuscaloosa,  AL 
35487-0350)  and  Robert  W.  Scharstein  (Electrical  Engineering  Department,  University  of 
Alabama,  Tuscaloosa,  AL  35487-0286) 

The  transmission  problem  of  time-harmonic  acoustic  scattering  by  the  velocity/density 
contrast  wedge  is  formulated,  through  symmetry  arguments  and  with  the  construction  of 
suitable  Green’s  functions,  as  a  pair  of  uncoupled  surface  integral  equations,  each  with 
one  unknown  function  dehned  on  a  single  wedge  face.  The  convergent  solution  to  the 
Fredholm  integral  equation  of  the  first  kind  is  e.Kpressed  as  a  Galerkin  series  of  Laguerre 
polynomials,  and  the  development  takes  account  of  the  distant  behavior  anticipated  from 
the  asymptotic  solution  to  the  related  Sommerfeld  half-space  problem.  The  required  inner 
products  for  the  Galerkin  projection  scheme,  which  are  integrals  of  products  of  the  weighted 
Laguerre  polynomials  and  the  pair  of  Green’s  functions  for  the  separate  homogeneous 
regions,  are  written  as  Taylor  series  coefficients  of  an  auxiliary  function.  Efficient  numerical 
implementation  of  the  physically-based  analysis  produces  an  accurate  and  intuitive  picture 
of  the  wave  interactions  with  this  canonical  scatterer.  [Work  supported  by  the  Naval 
Research  Laboratory  -  Stennis  Space  Center.] 
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Mellin  Transform  Solution  for  the  Static 
Line-Source  Excitation  of  a 
Dielectric  Wedge 

Robert  W.  Scharstein,  Member,  IEEE 


Abstract — An  integral  transform  analysis  of  the  static  scatter¬ 
ing  of  the  two-dimensional  potential  radiated  by  a  iine  source  in 
the  vicinity  of  a  penetrable  wedge  is  presented.  The  Mellin 
transform  is  used  to  derive  the  exact  static  solution  to  Laplace’s 
ti|uation  for  the  dielectric  wedge,  in  the  form  of  a  modal  series. 
The  important  dieiectric  edge  condition  behavior  is  explicitly 
contained  in  this  analytic  solution. 


I.  Introduction 

HE  static  solution  is  a  critical  component  of  any 
dynamic  scattering  problem,  especially  with  regard  to 
edge  conditions  and  behavior  in  source  regions.  A  qua¬ 
sistatic  philosophy  of  extracting  the  dominant  behavior  of 
a  solution  to  the  Helmholtz  equation  (V^  +  k~)tl/(r)  =  0 
from  the  static  (k  =  0)  solution  to  Laplace’s  equation 
V  T/f(f )  =  0  in  the  neighborhood  of  boundary  discontinu¬ 
ities  and  sources  is  well  known  and  successful  [l]-[4],  and 
is  naturally  advocated  by  Anderson  and  Solodukhov  [5] 
for  the  dielectric  wedge.  The  important  edge  behavior  of 
the  dynamic  fields  in  the  immediate  vicinity  of  the  apex 
r  =  0  is  obtained  from  the  static  analysis,  since  all  physical 
dimensions  are  scaled  by  the  wavelength,  and  /cr  -♦  0  is 
the  limit  of  interest.  The  numerical  results  of  Marx  [6] 
indicate  that  the  dynamic-held  behavior  at  the  apex  of  the 
dielectric  wedge  can  differ  from  the  static  solution.  How¬ 
ever,  continued  mathematical  analysis  of  the  wave  prob¬ 
lem  is  required  before  initiating  a  comprehensive  physical 
interpretation. 

This  paper  presents  the  complete  and  exact  solution  to 
the  static  dielectric  wedge  problem.  It  should  be  noted 
that  the  special  case  with  the  line  source  lying  on  the  x 
axis  (the  wedge  bisector)  does  appear  in  the  Russian  text 
[7].  Furthermore,  Smythe  [8]  presents  a  formal  treatment 
of  the  static  dielectric  wedge,  which  is  unfortunately  flawed 
with  divergent  integral  representations.  In  any  event,  the 
exact  static  solution  is  now  available  for  this  canonical 
geometry.  The  potential  problem  is  solved  using  the  Mellin 
transform,  which  is  itself  the  static  limit  of  the 
Kontorovich-Lebedev  transform  employed  by  Jones  [9] 
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for  the  scattering  of  time-harmonic  waves  from  impene¬ 
trable  wedges  having  Dirichlet  and  Neumann  boundary 
conditions.  The  static  modes  of  [5]  appear  as  the  residues 
in  the  Mellin  inversion  integral. 

II.  Statement  of  the  boundary- value  problem 

The  static  excitation  of  the  penetrable  wedge  is  due  to 
the  line  charge  of  unit  lineal  density,  located  at  the  source 
coordinates  (F,  <^')  of  Fig.  1.  The  permittivity  of  the 
wedge  of  angle  2  a  is  £,>  which  is  surrounded  by  a 
medium  with  permittivity  c,.  All  geometry  and  therefore 
all  resultant  flelds  are  of  inflnite  extent  and  invariant  in 
the  z  dimension,  which  restricts  the  phvsical  domain  to 

The  irrotational  electrostatic  fleld  is  uniquely  character¬ 
ized  by  the  scalar  potential  tpir),  which  is  a  solution  of 
Poisson’s  equation  subject  to  appropriate  boundary  condi¬ 
tions.  Denote  by  the  potential  fleld  in  the  external 
region  where  the  source  is 

,  1 

4,)= - 8(r  -  r')dU  -  4>') 

e,r 

{a  <  4>  <  2it  —  a),  (1) 

and  let  be  the  source-free  fleld  inside  the  wedge 

V^i/fjCr,  d»)  =  0  {-a<(i><a).  (2) 

Boundary  conditions  at  the  material  interfaces  4>  =  ±a 
for  this  scalar  potential  are  continuity  of  (from  continu¬ 
ity  of  the  tangential  electric  fleld)  and  continuity  of  the 
normal  electric  flux  density  e^^ipj/dn  (absence  of  any  free 
surface  charge). 

In  order  to  simplify  the  ensuing  analysis,  it  is  expedient 
to  decompose  the  desired  solution  i/r(x,  y)  for  the  bound- 
aiy-value  problem  of  Fig.  1  into  its  odd  4/°ix,  y)  and  even 
>P"(x,y)  components 

ipix,  ±y)  =  -[ip‘{x,y)  ±  il/°ix,y)]  (y  >  0)  (3) 

with  respect  to  the  x  axis  which  bisects  the  wedge.  Since 
ij/°{x,y)  is  an  odd  function  of  y,  it  vanishes  on  the  y  =  0 
plane  and  is  therefore  the  solution  in  the  upper  half-space 
y  ^  0  for  the  problem  having  a  soft  bisecting  plane 
(Dirichlet  boundary  condition).  This  is  equivalent  to  an 
out-of-phase  image  source.  Similarly,  (^'(.r,  y)  is  the  solu- 
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Fig.  1.  Dielectric  wedge  and  static  line  charge. 


while  the  discontinuity  in  normal  derivative 

—  4)')  -  <f>')  =  —Sir  -  r')  (10) 

d<p  oq>  Cl 

results  from  integrating  (1)  from  d>'  -  to  d>'  +  . 

The  Mellin  transform  [10]  is  used  to  transform  the 
radial  variable  r  to  a  complex  variable  s,  whereupon  the 
remaining  differential  equation  with  boundary  conditions 
in  is  solved  in  closed  form.  The  physical  solution 
ilfir,  4>)  is  recovered  by  careful  evaluation  of  the  inverse 
Mellin  transform.  The  Mellin  transform  of  f(r)  is 


F(s)  =  f  r*  ^ fir)  dr, 

J(\ 


and  has 


Fig.  2.  Odd  and  even  symmetry  components. 

tion  for  the  hard  bisecting  plane  (Neumann  boundary 
condition),  which  sustains  an  in-phase  image  source.  This 
symmetry  is  depicted  in  Fig.  2. 

III.  Mathematical  Analysis 
A.  Mellin  Transform  for  the  Case  of  Odd  Symmetry 

The  source  coordinate  4>  =  <t>'  divides  region  (1)  into 
two  source-free  regions,  resulting  in  a  total  of  three 
subregions  to  consider  for  the  half-space  above  the  per¬ 
fectly  soft  plane: 


fir)  =jrHFis)-,  r)  =  - — :  r  ‘  r  ’Fis)  ds  (12) 
Ztti-'c-/® 

as  its  inversion  formula.  If  r'‘~'fir)  is  absolutely  inte- 
grable  on  the  positive  real  axis  for  some  a  >  0,  then  the 
inversion  is  valid  for  c  >  a. 

The  Mellin  transform  of  the  product  of  and  Poisson's 
equation  (1) 

I  d  di-  \ 

=  -—Sir  -  r')Si<f>- <t>')  (13) 


(  4>2^r,  4>), 

0  £  <t>  ^  a 

is  the  simple  form 

a  <  <f>  ^  <f>' 

[>l't(r,<f>). 

4>'  <  <p  ^  rr. 

1 

The  Dirichlet  boundary  condition  on  the  soft  plane  is 

tl,2ir,0)  =  0  (5) 

i/rf(r,  it)=0  (6) 

and  at  the  material  interface 

iffzir,  oi)  =  «)  (7) 

<?  d 

Consistent  with  the  conventional  Green’s  function  ansatz, 
the  potential  is  continuous  everywhere  across  the  source 
plane  ()>=  4>' 

iAr(r,<^')  =  iAr(r,<^'),  (9) 

-sin(5a)  sin(sa)  cos  (sa) 

—  cos(sa)  — cos(sa)  sin(5a) 

Cl 

0  sinfjd'')  cos  is<i>') 

0  cosis4>')  -sin(5<^') 


This  fortunate  property  of  the  Mellin  transform  of  the  r 
dependence  in  the  two-dimensional  Laplacian  is  responsi¬ 
ble  for  its  successful  application  [7],  [10],  [11]  to  potential 
problems  in  wedge-shaped  regions.  In  region  (2)  where 
there  is  no  forcing  term,  this  procedure  gives 

+ -^W,(s,d>)  =  0.  (15) 


The  complex  variable  5  is  a  parameter  in  the  above  pair 
of  ordinary  differential  equations  in  4>,  with  solutions 

^2('t.  <(>)  =  ^(■t)sin  (sd»)  (0<d><a)  (16) 

■4'f(s,  <j>)  =  B(j)sin  is4>)  +  Cis)  cosis<f>) 

ia  <  <b  <  <t>')  (17) 
(s,  d>)  =  D(j)  sin  [s(d>  -  tt)]  i<t>' <  <f>  <  v)  (18) 

in  view  of  the  soft  boundary  conditions  (5)  and  (6).  Trans¬ 
formation  of  the  four  remaining  conditions  (7)-(10)  gives 
the  set  of  simultaneous  equations 

0  L  .  r  n  1 


sin  [j(7r  -  </)')] 
-cos[s(7r  -  <)>')] 
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to  be  solved  for  the  coefficient  functions  in  (16)-<18). 
Solution  of  these  yields  the  (soft)  Mellin  transforms 

1  -  r  r" 

- - sin[j(7r-  0')]sin(5<^)/A^(j) 

€,  S 

(0^<f>^a)  (20) 

r” 

^f(s,d>)  =  - sin[j(7r-  (f>')]{sin  (j<^) 

-t-r sin[jr((A  -  2ar)]}/Ao(j) 

(21) 

r'* 

^f(5,  (t>)  =  —  sin  [j(7r  -  <^)]{sin  (j<^') 

-f-r sin[5(<^'  -  2a)]}/Ag(s) 

<4>^  tt)  (22) 

with  denominator  function 

Agfj)  =  sin(s7r) -t- r  sinfjfTT  -  20;)],  (23) 

and  dielectric  contrast  parameter 

r  =  -^^.  (24) 

B.  Modifications  for  the  Case  of  Even  Symmetry 

In  the  case  of  a  hard  ground  plane,  the  Neumann 
boundary  conditions 


—  iA2(r,0)  -0 
d<p 

(25) 

d 

—  (A['(r,  77)  =  0 
d<p 

(26) 

replace  the  Dirichlet  boundary  conditions  (5)  and  (6)  of 
the  previous  section.  A  similar  application  of  the  Mellin 
transform  and  the  other  unchanged  boundary  conditions 
yields  the  (hard)  Mellin  transforms 

r  -  1  r'* 

'^As,4>)  = - cos [5(77-  (t>'')]cos(s<t>)/AXs) 

«i  ^ 

(Os0<a)  (27) 
r'* 

'Pr(5,  <f))  =  —  cos  [5(77  -  0')]{r cos  [5(0  -  2a)] 
e,j 

-cos(s<^)}/Aj(s),  {a  <  (t>  ^  <i>')  (28) 
r'* 

0)  =  - COS  [5(77  -  <^)]{r  cos  [sift'  -  2a)] 

-cos{s<l>'))/\,{s)  (4>'  <4>^  tt)  (29) 

where  in  this  case  the  denominator  function  is 

A, (5)  =  sin  (stt)  -  r  sin  [5(77  -  2a)].  (30) 

Note  that,  except  for  the  simple  scaling  by  c,  which 
persists  from  the  original  source  strength  chosen  in  (1), 


the  presence  of  two  different  dielectrics  is  entirely  ac¬ 
counted  for  by  r  in  all  of  the  above  transforms. 

C.  Inverse  Mellin  Transforms — Preliminaries 
The  zeros  of  the  denominator  functions  (23)  and  (30) 

A(i)  =  sin  (stt)  ±  r  sin  [5(77  -  2a)] 

are  central  to  the  Mellin  inversion  (12).  These  real,  simple 
zeros  can  be  computed  via  Muller’s  algorithm  [12]  for 
arbitrary  half-angle  a,  or  solved  analytically  as  the  roots 
of  a  trigonometric  pol5momial  when  a  is  a  rational  multi¬ 
ple  of  77.  This  procedure  is  demonstrated  for  the  particu¬ 
lar  case  a  =  77/3,  whereupon  the  variable  change  u  = 
sir/ 3  in  (31)  gives 

sin  (3«)  ±  r  sin  («)  =  0  (32) 

which  factors  into 

[3  ±  r  -  4  sin^  (u)]  sin  (u)  =  0.  (33) 

The  required  roots  are  now  explicitly  given  by 


(m  =  0,±  1,  ±  1,  ••• ),  (34) 

where  it  is  recalled  that  the  -t-  F  ( -  D  denotes  the  case  of 
soft  (hard)  symmetry.  The  effect  of  the  dielectric  material 
(e,  #  Cj  =»  r  #  0)  on  the  potential  above  both  symmetry 
planes  is  a  regular  displacement  of  the  integer  poles  for 
the  homogeneous  case  (e,  =  e,  =»  F  =  0).  The  index  n  in 
(34)  is  a  denumerable  ordering  of  these  poles. 

D.  Inverse  Mellin  Transform  for  the  Case  of  Odd  Symmetry 

As  r  -»  0  the  odd  potential  i/r"  -♦  0  and  the  dipole 
behavior  t/r®  ~  1/r  prevails  as  r  -*  tc  in  the  far  field.  A 
sufficient  Bromwich  contour  for  the  complex  integration 
(12)  is  therefore  guaranteed  for  the  choice  of  real  con¬ 
stant  0  <  c  <  1.  Complete  details  of  the  Mellin  inversion 
for  the  potential  <f))  of  (4)  inside  the  dielectric  sector 
are  provided,  whereupon  the  final  forms  for  ilfiir,  <f))  and 
ftiir,  <f>)  are  immediately  written  by  comparison. 

The  sin(jd>)  factor  in  'F,(j,  d>)  of  (20)  together  with  the 
transform  property  [10] 

{sin  (s(i>) Fis);  s]  =  -J^ifire'*)]  (35) 
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Fig.  3.  Bromwich  contour  in  the  complex  s  plane. 


renders 


f(r)  = 


sin  [s(7r  -  <f>')](r'/ry 
s{sin  (sir)  +  f  sin  [rCrr  -  2a)]} 


the  desired  function.  The  integration  path  in  (36)  has 
been  pushed  flush  against  the  imaginary  axis,  and  the 
principal  value  notation  invoked  to  properly  account  for 
the  pole  at  the  origin.  For  r  >  r',  closure  at  infinity  in  the 
right  half-plane  gives  a  convenient  contour  on  which  to 
apply  Cauchy’s  integral  theorem,  and  is  shown  in  Fig.  3. 
Let  C„  be  a  circle  of  vanishing  radius  p  -»  0,  centered  on 
the  pole  s„  with  residue 

lim  —6  G(s)ds  - -  ,  - ,  (37) 

p->0  2vi~C„ 

where  the  derivative  of  the  soft  denominator  function  is 
cos(j„7r)  -I-  (it  -  2a)r cos  [j„(7r  -  2a)]. 


The  integral  around  the  origin 

hm  —6  G(s)  ds  =  -  ■  (39) 

p-.o27ri~Co  A'„(0) 

follows  from  the  limit  s„  -*  0.  Let  positive  n  =  1,2,  be 
the  indices  of  the  poles  in  the  right  half-plane,  and  let 
negative  n=  -1, -2,  identify  the  poles  in  the  left 
half-plane.  The  odd  symmetry  of  the  function  A^fs)  pro¬ 
vides 

=  and  A'„(-r,)  =  A',(r„).  (40) 

The  contribution  from  the  integral  around  the  infinite 
semici'cle  R  of  Fig.  3  is  zero  for  r  >  r',  while  closure  of 
the  contour  in  the  left  half-plane  is  appropriate  when 
r  <  r' .  Cauchy’s  integral  theorem  now  yields 


“  sin  [5„(7r  -  <f>')]  ^ 


(r^r),  (41) 

which  together  with  (35)  gives 

r  -  1 

ip2(r,  <t>)  =  - 

^1 

“  sin[j„(7r-<^>')]sin(r>) 

■  1- - ' 

n-I 

ir^r')  for  the  sector  0  <  s  a  (42) 

as  the  inverse  Mellin  transform  of  (20).  Since  all  three  of 
the  transforms  (20)-(22)  are  of  the  same  general  form, 
the  remaining  two  field  expressions  are  evidently 

-1  *  sin  [5,(77  -  <f>')] 

^-ir,cf>)  =  —  I  -V.-  -■ 

Cl  „-i  2,A„(s,) 

•{sin(5,(^)  +  r  sin  [5„(d>  -  2a)]} 

(r  $  r')  for  the  seaor  a  <  4»  <  <b'  (43) 


- 1  *  sin  [5,(77  -  d))] 

Z  — TFcT)— ' 

'1  n-1 

•{sin  (s„t^>')  +  r  sin  [5,(d>'  -  2a)]} 

(r  $  r')  for  the  sector  d>’  <  (i>  <  77.  (44) 

Note  the  reciprocity  in  the  last  two  expressions  for  d>  ♦-*  </>'. 
The  above  reduce  to  the  correct  potential  due  to  a  unit 
line  charge  located  at  </>'  =  77/2  above  a  soft  ground 
plane  in  the  special  case  of  no  dielectric  wedge  (c,  =  c,). 

£.  Inverse  Mellin  Transform  for  the  Case  of  Even  Symmetry 

The  monopole  potential  In  r  due  to  the  original  line 
source  and  the  in-phase  image  is  the  dominant  feature  of 
ij/*  in  the  far  field  as  r  -♦  tc.  Therefore,  the  condition  of 
integrability  following  the  transform  pair  of  (11)  and  (12) 
is  not  satisfied,  and  the  Mellin  inversion  formula  cannot 
be  directly  applied  to  the  functions  in  (27)-(29).  Mathe¬ 
matically,  operating  with  the  rf)  derivative  of  the  potential 
temporarily  removes  this  troublesome  logarithmic  varia¬ 
tion.  which  is  then  restored  following  the  convergent 
Mellin  inversion. 

The  (t>  derivative  of  (27) 

,  1  -  r  cos  [5(77  -  di')]  sin  (5(f>)  . 

—%{s,cl>)  = - — - - r'’  (43) 

d<f>  e,  A, (5) 

is  transformed  in  the  manner  of  the  previous  section  to 

<?  r  -  1 

—  ^2^r,<j>)  = - 

dtp  e, 

*  cos[5,(77  -  (^')]sin(5„<f>)  ,, 

■  ^ 
n-1 

{r  $  r'),  (46) 

where  the  s„  are  now  the  hard  poles  in  (34)  and  -T 
replaces  F  in  (38)  for  this  hard  denominator  function  (31)- 
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The  antiderivative  with  respect  to  yields 

1-r  “  cos[j,(7r- d>')]cos(j>) 

* -  2-  - ^77—:;; - (r/r'r  ■ 

«i  ,-i 

(r  $  r')  for  the  sector  0  ^  ^  a.  (47) 

where  the  boundedness  at  r  =  0  and  the  known  behavior 
at  r  -*  ao  specify 

(  0,  for  r  <  r' 

«'?('•)  =  I  _1  In  (r/r’),  for  r  >  r' 

Tre, 


as  the  (^independent  term  in  this  solution  of  Laplace’s 
equation.  Similarly,  the  spatial  potentials  from  the  related 
transforms  (28)  and  (29)  are  thus 


-1  ”  cos  [s,(7r  -  <^')] 

1^;  (r,  (f»)  =  ^F^(r) - £  - -  X - ('■A')  ’ 


-l  n-l 


•{r cos[j„(</»  -  2a)]  -  cos(s„<f»)} 

(r  $  /)  for  the  sector  a  <  <f>  (f)'  (49) 


and 


-1  ”  cos[s„(7r  -  </))] 

=  w^(r)  -  —  Z — tTlTI — ' 

'1  n-l 

•{rcos[j„(<f>'  -2a)]  -  cosis„4>’)} 
ir  $  r')  for  the  sector  <f»'  <  ^  ir.  (50) 

These  results  are  also  verified  for  the  special  case  of  a 
unit  line  charge  located  at  <f>'  =  -ir/2  in  a  homogeneous 
half-space  (^2  =  ^i)  above  the  hard  symmetry  plane. 


IV.  Results  for  the  Complete  Static  Solution 

Contours  of  constant  potential  are  illustrated  in  Fig.  4 
for  the  case  of  wedge  half-angle  a  =  tt/S.  Even  and  odd 
potentials  are  computed  for  the  separate  boundary-value 
problems  and  combined  according  to  (3).  The  line  source 
is  located  a  unit  distance  (r'  =  1)  from  the  wedge  apex  at 
the  {x,  y)  coordinate  origin.  The  domain  of  Fig.  4  is  ail 
internal  to  the  unit  circle.  Note  that  in  this  static  problem 
with  material  boundaries  at  <f>  =  ±a,  the  source  distance 
r'  is  the  single  length  scale.  As  62/^1  large,  the 
dielectric  wedge  becomes  essentially  an  equipotential  re¬ 
gion,  which  is  expected  since  an  increasingly  dense  dielec¬ 
tric  (62  -^  “)  behaves  electrically  like  a  perfect  conductor. 
As  r  -►  0,  the  dominant  behavior  of  the  scalar  potential  is 
given  by  the  first  terms  of  the  modal  series  (42)-(44)  and 
(47)-(50),  as  pointed  out  by  Anderson  and  Solodukhov  [5] 
and  tabulated  in  the  recent  text  by  Van  Bladel  [13]. 


V.  Conclusions 

The  exact,  modal  series  solution  of  Laplace’s  equation 
for  the  static  field  inside  and  outside  of  the  dielectric 
wedge  is  obtained  via  the  Mellin  transform,  and  is  readily 
evaluated.  This  analytic  solution  is  a  useful  benchmark 
with  which  to  check  any  numerical  and  approximate 
methods  that  are  applied  to  the  static  problem. 


Fig.  4.  Contours  of  constant  potential  for  the  dielectric  wedge.  Case: 
r'  -  1.  (6'  -  77/2.  a  -  77/3,  e,  =*  1,  e,  =  10. 
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